Exotic d-wave superconductivity in strongly hole-doped K rr Bai_ a; Fe2As2 
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We investigate the superconducting phase in the K a: Bai_ a; Fe2As2 122 compounds from moderate 
to strong hole-doping regimes. Using functional renormalization group, we show that while the 
system develops a nodeless anisotropic s ± order parameter in the moderately doped regime, gapping 
out the electron pockets at strong hole doping drives the system into a nodal (cos k x +cos k y )(cos k x — 
cosky) d-wave superconducting state. This is in accordance with recent experimental evidence from 
measurements on KFe2As2 which observe a nodal order parameter in the extreme doping regime. 
The magnetic instability is strongly suppressed. 

PACS numbers: 74.20.Mn, 74.20.Rp, 74.25.Jb, 74.72. Jb 
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The most elementary questions in the field of iron- 
based superconductors, such as the symmetry of the or- 
der parameter in the superconducting (SC) state, are 
still under vivid debate. The complexities involve an 
intricate band structure, a diversity of different mate- 
rial compounds which exhibit sometimes contradictory 
behavior, and the proximity of various symmetry-broken 
phases. Due to best single-crystal quality, the most stud- 
ied pnictidc compounds belong to the 122 family such as 
BaFe2 As2- Their crystal structure is tetragonal I4/mmm, 
where the Fe and As atoms arrange into layers; the intra- 
layer hybridization is dominant but, unlike other pnictidc 
compounds such as the 1111 family, the inter- layer hy- 
bridization is also important. Soon after their discov- 
ery (Hi the 122 compounds have been synthesized not 
only with Ba as a substituent between the FeAs layers, 
but also with K, Cs, and Sr. The SC transition temper- 
atures achieved were up to 37 K Q. 

The current theoretical opinion on the SC order pa- 
rameter has converged on a nodeless order parameter 
that changes sign between the electron (e) pockets and 
hole (h) pockets. This order parameter comes out of 
both the strong and the weak-coupling pictures of the 
iron-based superconductors 0-0], and owes its origin to 
the pnictide Fermi surface (FS) topology of h pockets at 
the T and e pockets at the X (7r,0)/(0,7r) point of the 
unfolded Brillouin zone. The dominant scattering con- 
tributions originate from h pocket scattering at V to e 
pockets at X, yielding the SC order parameter for 
the doped case and the collinear antifcrromagnctic phase 
in the undoped case. Detailed nesting properties of the 
pockets, the multi-orbital character of the FS, and the 
presence or absence of a third h pocket at M (ir, it) in 
the unfolded Brillouin zone complicate this picture. For 
the 1111 compounds, it was shown that the absence of 
the M h pocket (whose Fermi level can be significantly 
tuned by the pnictogen height through replacing As by 
P) can modify the SC order parameter anisotropy from 



a nodeless to a nodal s phase, which gives the correct 
material trend for As-P substitution in other pnictidc 
families (zHsl- With small exceptions, the anisotropic ex- 
tended s-wave scenario (and its extension to the nodal 
s^) was consistent with experimental findings for most 
of the pnictide compounds Tol-Hij . 

It was realized at a very early stage that electron 
and hole doping can have qualitatively different effects 
in the pnictides [lj|. Hole doping should increase the 
propensity to a nodeless (s ± ) SC phase. The qualita- 
tive picture applies to both the 122 and the 1111 com- 
pounds: as the Fermi level is lowered, the M h pocket 
becomes more relevant and the M -f> X scattering adds 
to the (ir, 0)/(0, 7t) scattering from T to X. As such 
the anisotropy-driving scattering such as inter-electron 
pocket scattering becomes less relevant and it yields a 
nodeless, less anisotropic and more stable |3.ll7j. This 
picture is qualitatively confirmed by experiments. While 
thermoelectric, transport and specific heat measurement 
have been performed for K 2 ,Bai_ a; Fe2As2 from x = to 
the strongly hole-doped case x = 1 H|-20j], more de- 
tailed studies have previously focused on the optimally 
doped case x = 0.4 with T c = 37 K, where all measure- 
ments such as penetration depth and thermal conductiv- 
ity find indication for a moderately anisotropic nodeless 
-H. ARPES on doped BaFe 2 As 2 , likewise, 
~-2t 



gap 



□J, 21 



finds a nodeless SC gap [2 



The experimental findings for the SC phase in KFe2As2 
came as a sur pris e. Thermal conductivity [27, 28 1, pene- 
tration depth [29( | , and NMR 3^, 31 1 provide clear indica- 
tion for nodal SC. The critical temperature for KFe2As2 
is at T c ~ 3 K, one order of magnitude less than the 
optimally doped samples. ARPES measurements [32[ 
show that the e pockets have nearly vanished, while the 
h pockets at the folded T point are large and have a lin- 
ear dimension close to ir/a. In this Letter, we provide a 
detailed picture of how the SC phase evolves under hole 
doping in K I ,Bai_ 2 ,Fe2As2, and find that the nodal phase 
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FIG. 1. (color online) (a) Schematic plot of the unfolded 
FS. (b) fc 2 = slice of the 122 band structure given in [3^ |. 
The dominant orbital weights along the e and h pockets are 
highlighted. The patches along pockets are enumerated coun- 
terclockwise, starting at each pocket with the patch indicated 
by a dot. The number of patches with electron pockets is 80, 
and without it decreases to 48. 



observed for x = 1 is of (extended) d-wave type. We use 
functional renormalization group (FRG) to investigate 
how the SC form factor evolves under doping from node- 
less anisotropic in the moderately hole-doped regime 
to ci-wave in the strongly hole-doped regime, where the e 
pockets are assumed to be gapped out. The d-wave SC 
minimizes the on-pockct hole interaction energy. We find 
the critical divergence scale to be an order of magnitude 
lower than for the optimally doped scenario, which is 
consistent with experimental evidence. 

We focus on studying K a; Bai_ :r Fe2As2 starting at the 
optimally doped case around x = 0.4 and increasing the 
hole doping up to KFe2As2- We use an effective 5-band 
tight-binding model developed by Graser et al. [33| to 
describe the band structure of the 122-type iron-based 
superconductors [see Fig. Q]: 
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Here c's denote electron annihilation operators, a, b the 
five Fe <i-orbitals, K the band matrix, and s the spin in- 
dex. As seen in Fig. Q] and Fig. [2] for moderate hole dop- 
ing, the conventional five pocket scenario with e pockets 
at X(tt, 0) and M(7r,7r) emerges. For larger hole dop- 
ing, the e pockets vanish and only small disconnected 
lobe features are found around X [Fig. dt] . The kinetic 
model reduces to the effectively three h pocket scenario 
shown in Fig.^. Other details of the 122 band structure 
are currently still under debate, with unresolved ques- 
tions about the FS topology at the Z point in the three 
dimensional Brillouin zone and the importance of inte- 
grating over the full range along k z . However, as many 
of these details mostly affect the e pocket anisotropics, 
they are irrelevant for our proposed SC mechanism: as 
we always consider a rather largely hole-doped regime, 
the e pockets can be assumed relatively small - even dis- 
appearing in the most interesting case, i. e. that of full 



hole doping. We, therefore, particularize to the k z = 
cut of ([1]) in the following, and also omit the lobe features 
at large hole doping within the RG calculations. To test 
our assumption of the irrelevance of the k z dispersion to 
our results, we have made several other cuts at different 
k z and confirmed that our results do not change qual- 
itatively. We cannot ultimately exclude that the lobes 
may influence the system due to the fact that our Bril- 
louin zone patching scheme is not fully adequate for such 
small Fermi surface features. Still, within our formalism, 
we find that the lobes are negligible in the RG flow. 

A schematic picture of the FS topology is given in 
Fig. [TJl. The h pockets at T mainly have d xz and d yz 
orbital content while the h pocket at M consists of d xy 
orbital weight. When present, the e pockets consist of 
d xz and d yz orbital weight. Exceeding a certain size, the 
front tip along T — X also has an important d xy weight 
on the e pockets. We use the conventional onsite orbital 
model for the interactions, i. e. 
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where n^ as denote density operators at site i of spin s 
in orbital a. We consider intra- and inter-orbital inter- 
actions Ui and U 2 as well as Hund's coupling Jh and 
pair hopping J pa i r . We choose the values of the interac- 
tion parameters close to the ones obtained by constrained 
RPA ab initio calculations [34J : U\ > U2 > Jh ~ ■/pain 
and set Ui = S.QeV, U 2 = 2.0eV, J H = J pa ir = O.GeV. 
While there are variations of these parameters for differ- 
ent classes of pnictides, the values of the parameters are 
all in the same range, and we have confirmed that vari- 
ations of 20-30 % of the interaction parameters do not 
change the picture qualitatively. As a tendency, a com- 
parably large absolute value of U\ needs to be kept to 
trigger the SC instability, where increasing U2 also helps 
to increase the critical cutoff scale and, thus, T c . 
Using multi-band FRG jllMS. 

we study the evo- 
lution of the renormalized interaction described by the 
4-point function (4PF) under integrating out high en- 
ergy fermionic modes: V A (k 1 ;k 2 ; fc 3 ; k 4 )cl 4S cl 3S c k2S c ki - s7 
where the flow parameter is the IR cutoff A approaching 
the FS. ki to k 4 are the incoming and outgoing mo- 
menta. The starting conditions are given by the bare 
initial interactions for the 4PF with the bandwidth serv- 
ing as an UV cutoff. The diverging channels of the 4PF 
under the flow to the FS signal the nature of the instabil- 
ity, and the corresponding A c serves as an upper bound 
for the transition temperature T c . The Cooper channel 
of the 4PF provides the different SC form factors - the 
dominant order parameter having the largest eigenvalue 
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FIG. 2. (color online) Representative scenarios of the FS (unfolded BZ) and instability eigenvalue flows for chemical potential 
and electron concentration per iron fi = — .12, n e \ = 5.913 (a), /x = —.22, n e i = 5.663 (b), and /i = —.32, n e i — 5.346 in (c). The 
hole doping of our model calculation in (c), while exceeding the experimental setup n a \ — 5.5, best matches the FS profile from 
ARPES [33I • The dominant and subdominant scatterings in the Cooper channel are highlighted in (al)-(cl) by full and dashed 
arrows. The color contours along the FS label the dominant orbital weights [inset (al)]. The leading eigenvalue flow of the 
ordering channel for different Fermi instabilities [charge density wave (CDW), Pomeranchuk (PI), spin density wave (SDW) 
and superconductivity (SC)] are plotted in (a2)-(c2) versus the IR cutoff FRG flow parameter A. For (a) and (b) we find s ± 
as the leading Fermi instability. For (c) we observe a leading d-wave instability. 



HSi^. In Fig.[U the leading eigenvalues for different 
FS instabilities are plotted against A for different fillings 
between moderately hole-doped from the left to strongly 
hole-doped to the right. We find that for all scenarios 
the leading instability is in the Cooper channel. 

For the moderately doped case, the e pockets arc of 
similar size as the h pockets. Fig. [2] (al) shows the FS 
structure as well as the dominant (full line) and subdom- 
inant scattering (dashed arrow) processes in the Cooper 
channel. The two major components are given byTf)-! 
as well asMol scatterings. They are particularly im- 
portant for the front tips of the e pockets since these 
parts can scatter to M via dominant U\ interaction due 
to identical orbital content. The SDW fluctuations are 
strong, signaling the proximity to the leading magnetic 
instability scenario of the undoped model [Fig. [2] (a2)]. 

For the intermediate regime, between moderate and 
strong hole doping, the e pockets are already very small 
[Fig. [2] (bl)]. The nesting to the h pocket is absent, 
and the SDW fluctuations are strongly reduced. In addi- 
tion, the SDW fluctuations become less concentrated in 
the (71", 0)/(0, 7r) or (71", 7r) channel, and spread into vari- 
ous incommensurate sectors [38| . The d xy orbital weight 
on the e pocket is reduced and the M <H- X scattering 
becomes subdominant. The main Cooper channel scat- 
tering is along r o X. As a consequence, s± is still the 
leading instability, where the form factor and its decom- 
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FIG. 3. (color online) (a) Form factor of the leading SC insta- 
bility of scenario Fig. [2b, plotted versus the patching index 
of the Fermi pockets according to Fig. [T^i. (b) Eigenvalues of 
the orbital decomposition of the SC form factor in (a). The 
ratio of the values label the relative importance of the orbital 
scattering channel V(a, a 
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position into orbital scattering contributions are shown 
in Fig. [3] the largest gap is found for the inner h pocket 
at r, followed by the outer h pocket and the h pocket 
at M, where the e pockets show anisotropic gaps. The 
orbital decomposition confirms the previous discussion 
of the the dominant scattering contribution, in that the 
largest weight resides at intra and inter-orbital scattering 
of the d xz and d yz orbital. However, we already observe 
that, due to the lack of SDW fluctuations supporting the 
SC, the critical divergence scale is decreased [Fig.[2](a2)- 
(c2)]. In particular, while still subdominant, we can al- 
ready see the d-wave evolving as the second-highest insta- 
bility eigenvalue in the Cooper channel. When e pockets 
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FIG. 4. (color online) (a) Form factor of the leading SC in- 
stability of the strongly hole-doped scenario Fig. [2b, plotted 
versus the patching index of the FS according to Fig. (b) 
Eigenvalues of the orbital decomposition of the SC form fac- 
tor in (a). We find a extended d-wave form factor. The nodal 
points are given along the main diagonal of the Brillouin zone, 
i.e. the form factor of the inner h pocket, for example, crosses 
zero between the patches (2, 3), (6, 7), (10, 11), and (14, 15). 



are still present, the form factor (not shown here) closely 
resembles the extended d-wave type involving h pockets 
and e pockets [3t| . 

At strong hole doping, the e pockets are absent, and 
the h pockets are very large. The flow in Fig. [2fc shows no 
instability up to small cutoff scales of A where we find a 
leading instability in the Cooper channel. Its form factor 
and orbital scattering decomposition is shown in Fig. [4] 
We observe an extended d-wave instability on the three 
h pockets, with nodes located along the main diagonals 
in the Brillouin zones [as seen by comparing the patch 
numbers of the crossings in Fig. HJt and the patch- 
ing enumeration defined in Fig. [T]. A harmonic analy- 
sis of the order parameter yields a large contribution of 
cos(2fca;) — cos(2fc y ) type and a subdominant cos(A; x ) — 
cos(fcy) component, i.e. the form factor is most accu- 
rately characterized as (cosfc^ + cos k y )(cos k x — cosk y ). 
The dominant scattering is intra-pocket scattering on the 
large M h pocket, followed by inter-orbital d xy to d xz ^ yz 
scattering between MoT. While the magnetic fluctu- 
ations are generally weak in this regime, the dominant 
contribution is now given by (7r,7r) SDW fluctuations 
as opposed to (7r, 0)/(7t, 0) for smaller hole doping. For 
strong hole doping, the h pocket at M is large enough to 
induce higher harmonic ci-wave SC through intra-pocket 
scattering between the d xy orbitals as confirmed by the 
large value of d xy — d xy pairing [Fig. 0J}]. Via scatter- 
ing to the other pockets, the SC instability is likewise 
induced there, however, with smaller amplitude than for 
the M pocket [Fig. Hk] . As opposed to conventional first 
harmonic cZ-wave, there is no sign change between the 
extended d-wave form factor on the M pocket and the T 
pocket according to cos^fc^) — cos(2fc y ) [Fig.0]. This pic- 
ture of a fc-space proximity effect from the M pocket to 
the F pockets is substantiated by our checks with calcula- 
tions involving the M pocket only, where we see a similar 
evolution of an SC instability (the divergence is lower, as 
the inter-orbital scatterings in the 3-pocket scenario help 
to renormalize the repulsive Coulomb interactions) . This 



matches the orbital decomposition of the SC form factor 
in Fig. |4Jd, showing dominant intra-orbital scattering of 
the d xy orbital. 

As apparent from the ARPES data, the nodal char- 
acter of the SC phase in KFe2As2 cannot originate from 
possible nodes on the e pockets (which are gapped out 
at these doping levels) but must be due to nodes on the 
h pockets. It is then clear that the order parameter can- 
not be as it does not tend to allow for an anisotropy 
that would drive the h pockets nodal. The d-wave insta- 
bility which we find for the strongly hole-doped regime 
provides an explanation for the general experimental evi- 
dence, while the detailed gap structure certainly deserves 
further investigation [39| . Electron- phonon coupling may 
change the picture slightly quantitatively, but not qual- 
itatively, as the nodal features tentatively linked to the 
d-wave symmetry are unambiguously observed in experi- 
ment. Pnictogen height variations as a function of doping 
may change the precise value of T c , and would be impor- 
tant to be studied in general from first principles. Finally, 
it would be interesting to further analyze how the system 
evolves from the s-wave SC phase to the d-wave SC phase 
as a function of doping. 
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